Ion-ion correlation and charge reversal at titrating solid interfaces 



OO 
O 

o 

(N 

o 

Q 



o 



Christophe Labbez 1 , Bo Jonsson 2 , Michal Skarba 3 and Michal Borkovec 3 
1 Institut Carnot de Bourgogne, UMR 5209 CNRS, 
Universite de Bourgogne, F-21078 Dijon, France 
2 Department of Theoretical Chemistry, Chemical Center, 
Lund University, FOB 124, S-221 00 Lund, Sweden 
A Department of Inorgarnic, Analytical and Applied Chemistry, Geneva University, 1211 Geneva J h Switzerland 

(Dated: December 1, 2008) 

Confronting grand canonical titration Monte Carlo simulations (MC) with recently published 
titration and charge reversal (CR) experiments on silica surfaces by Dove et al. [lj and van der 
Heyden et al. we show that ion-ion correlations quantitatively explain why divalent counterions 
strongly promote surface charge which, in turn, eventually causes a charge reversal (CR). Titration 
and CR results from simulations and experiments are in excellent agreement without any fitting 
parameters. This is the first unambiguous evidence that ion-ion correlations are instrumental in 
the creation of highly charged surfaces and responsible for their CR. Finally, we show that charge 
correlations result in "anomalous" charge regulation in strongly coupled conditions in qualitative 
desagreement with its classical treatment. 

PACS numbers: 82.30.Fi, 82.20.wt, 82.65.+r, 82.70.Dd, 05.70.Np 
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A key parameter for understanding electrostatic inter- 
actions in charged suspensions is the surface charge den- 
sity, a [H, IH ■ This parameter may control phase separa- 
tion, phase transition [H| , promote dissolution Q , control 
particle growth Ml, colloidal stability 0, 0| and mem- 
brane selectivity [10( and play a significant role in catal- 
ysis and fluid transport through ionic channels 
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Typically a surface becomes charged through a titration 
process where surface groups ionize, e.g. the titration of 
silanol groups -Si-OH -> -Si-0~ + H+. 

Ion-ion correlation in solutions of charged colloids has 
attracted much interest due to its ability to generate net 
attractive interaction between equally charged colloids. 
This is, however, only one aspect of ion-ion correlation 
and an equally important aspect is its facilitation of sur- 
face titration processes. An example of vital importance 
is the surface titration of calcium silicate hydrates in ce- 
ment paste, which reach a very high surface charge den- 
sity due to ion-ion correlation among the calcium counte- 
rions. The high surface ch arg e density is then responsible 
for the setting of cement [HI, 0, EH . 

An additional effect, due to ion-ion correlation, is when 
a charged surface immersed in a multivalent electrolyte 
more than compensates its charge resulting in a charge 
reversal (CR) |16lll7|. That is, the surface attracts coun- 



terions in excess of its own nominal charge and the ap- 
parent charge seen a few Angstrom from the surface ap- 
pears to be of opposite sign to the bare surface. CR 
is experimentally characterized by several methods, e.g. 
streaming current [2| and can induce repulsion between 
oppositely charged surfaces ID, 19j . Despite a vast num- 
ber of studies the origin of charge reversal is still con- 
troversial as highlighted in a recent review article [2fj| . 
where Lyklema points out the lack of discrimination be- 



tween the chemical (specific chemical adsorption) and 
physical (Coulombic interactions including charge corre- 
lations) origin of CR. 

The charging of interfaces (mineral or organic) has 
traditionally been described by the surface complexa- 
tion model (SC) [H], [SI, where titration of sur- 
face groups is described by an ensemble of mass balance 
equilibria with associated equilibrium constants. The 
SC is based on the Poisson-Boltzmann (PB) equation 
augmented with a so-called Stern layer. The PB equa- 
tion, however, neglects ion-ion correlation and e.g. fails 
to predict the attraction between equally charged sur- 
faces 13, 25| as well as the repulsion between oppositely 
charged ones [lj|. Despite its widespread use the limi- 
tations of the SC model in highly coupled systems have 
never been examined. Moreover, although multivalent 
ions are ubiquitous in many systems, there are only few 
experimental studies with titration of metal oxide parti- 
cles in the presence of such ions |l|, l2fj, l27[ . 

Recently a grand canonical titration simulation 
method was developed [28| and applied to a calcium sil- 
icate hydrate (C-S-H) solution [l3(. An excellent agree- 
ment was found between experiment and simulation for 
the titration process as well as for the CR. Unfortunately, 
due to the solubility range of C-S-H, the full titration 
curve is not accessible experimentally. 

In the present letter, we scrutinize the prediction from 
our MC simulations for the charging process and CR by 
comparing to experimental data for silica particles dis- 
persed in sodium and calcium chloride solutions [lj. In 
particular, we focus on the contribution of the correla- 
tions (ion- ion and ion-site). We also present a simu- 
lated charge regulation of two parallel silica surfaces in 
a Ca 2+ salt solution when decreasing their separation. 
The experimental investigation of surface titration of sil- 
ica fume particles by Dove et al. [1] and CR by multi- 
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valent counterions from a streaming current analysis in 
silica nanochannels by van der Heyden et al. will serve 
as references. 

The simulations are based on the Primitive Model with 
the solvent as a structureless medium characterized by 
its relative dielectric permittivity e r . All charged species 
are treated as charged hard spheres and the interaction, 
between two charges i and j separated a distance r is, 

u(r) = p^- if r>d hc (1) 

otherwise u(r) = oo. Z{ is the ion valency, e the elemen- 
tary charge, eo the dielectric permittivity of vacuum and 
dhc the hard sphere diameter of an ion always equal to 
4 A. The particles are modeled as infinite planar walls 
with explicit titratable sites with the intrinsic dissocia- 
tion constant Kq distributed on a square lattice. The 
titration process can be described as, 

MOH ^ MO- + H+, K Q = (2) 

a-MOH 

where the a's denote the activities of the species. The 
model is solved using a Grand Canonical Monte Carlo 
method (29[, i.e. at constant pH and chemical potential 
of the ions. In addition to ordinary MC moves the titrat- 
able sites are allowed to change their charge status. We 
imagine the deprotonation of a surface site as a two-step 
process: the release of a proton from the surface followed 
by an exchange of an ion pair (H + , B~) with the bulk. 
The corresponding Boltzmann factor for the trial energy 
can be expressed as, 

Nr 

exp(-/3AU) = -^exp[-/?0u B +AE7 ei )+mlO(pH-p# o )] 

(3) 

where fj, represents the chemical potential of a particular 
ion, V is the volume, Nb the number of ions B~ and 
AU el the change in electrostatic energy. An analogous 
expression holds for protonation, for more details see [l3j . 
To establish the model parameters for a silica/solution 
interface, we turn our attention to the case of silica in 
1:1 salt solution. The following set of parameters: pKo 
= 7.7, p s = 4.8/nm 2 , the surface site density, and dj s = 
3.5 A, the minimum separation between an ion and a 
surface site, have been found to give a perfect description 
of the charging process, see FigfTJ The surface charge 
density, a increases, in absolute value, with pH, but is 
always smaller than in the ideal case due to the strong 
electrostatic repulsion between the deprotonated silanol 
groups. A change in the salt concentration modulates 
the screening of the electrostatic interactions which, in 
turn, results in a change in a. 

When the sodium counterions are replaced by divalent 
calcium ions, the simulations predict a strong increase 
in the surface charge density of silica in perfect agree- 
ment with the titration data of Dove et al. Ill - see 




FIG. 1: Simulated (MC) and experimental (symbols) [l| sur- 
face charge density versus pH for silica particles dispersed in 
a sodium (1:1) salt solution at different concentrations. The 
ideal titration curve (ideal) is given as a reference 
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FIG. 2: Simulated (MC) and experimental (symbols) [l|] sur- 
face charge density versus pH for silica particles dispersed in 
a calcium (2:1) and sodium (1:1) salt solution at various con- 
centrations. The corresponding simulated titration curves for 
particles having implicit sites (MC imp.) are also given for 
comparison 

Fig|2] A similar behavior has been observed from titra- 
tion experiments on titania particles [27j in multivalent 
electrolytes and, indirectly, from proton release measure- 
ments following adsorption of multivalent ions (or poly- 
electrolytes) [3fJ emphasizing the universality of the phe- 
nomenon. The classical SC approach fails to predict the 
effect, unless a complexation reaction of calcium with 
the silanol surface groups is invoked. Here, we unam- 
biguously show that the charge correlations (ion-ion and 
ion-site), without any additional fitting parameters, are 
enough to explain the increased silanol group ionization 
in the presence of Ca 2+ ions. Thus, the complexation 
constants often invoked in the SC has no physical basis. 
In an attempt to distinguish the effect of ion-ion from 
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ion-site correlations a calculation was performed using a 
smeared out surface site density - see Figl^l A significant 
drop in |<j| is observed, which becomes more pronounced 
when the counterion valence is high, although the qual- 
itative behavior of a upon increasing the counterion va- 
lency remains the same. Thus, the discrete nature of the 
surface groups have a quantitative but not qualitative 
effect and we can conclude that ion-ion correlations are 
more important than ion-site correlations, in contradic- 
tion with recent approximate studies 
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FIG. 3: Comparison between simulated (line) and experimen- 
tal (points) 2J streaming currents of a silica channel as a func- 
tion of electrolyte concentration of a CaCb and MgCb salt 
at pH 7.5. In the simulations, the site density is 4.8 nm -2 , 
pK = 7.7, d ls (Ca 2+ ) = 3.5 A, d ls (Mg 2+ ) = 3.1 A. 

We can also compare our simulation data to the stream- 
ing current, I s tr, experiments through silica channels per- 
formed by van der Heyden et al. [2j in presence of di- and 
tri-valent counterions at pH 7.5. I s t r is computed from 
the simulated local charge density, p(x). For a channel of 
height h and width w with h >> w, I s t r can be defined 
as. 



2w 



h/2 



p(X)v(X)dX with X = x 



(4) 



where x is the distance normal to the silica surface 
and v{x) is the local fluid velocity, approximated as a 
Poiseuille flow with boundary conditions: v(X = 0) = 
and v'(X — h/2) — 0. The zero velocity is arbitrar- 
ily set at x = 3/2e?/j C where by definition we also find 
the electrokinetic potential. The simulations predict for 
both Mg 2+ and Ca 2+ a monotonic decrease in I str upon 
increasing salt concentration as also found experimen- 
tally - see FigfJJ A sign reversal in the streaming po- 
tential is predicted at ccr ~ 310 mM for Ca 2+ and 
ccr « 280 mM for Mg 2+ in good agreement numbers 
observed experimentally Q. What is more, the mag- 
nesium value is predicted to be below that of calcium 
in agreement with experiments 0] due to a higher sil- 
ica charge. Simulations were also performed for trivalent 



cations, where CR is predicted at sub-millimolar concen- 
trations, ccr ~ 400/iM, although times higher than what 
was observed with cobalt(lll)sepulchrate [2|. 
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FIG. 4: The degree of ionization, a, as a function of surface 
separation at various pH values. The site density is 4.8 nm~ 2 
and pKo = 7.7 and the surfaces are in contact with a 0.2 mM 
CaX2 solution. 

We now investigate the surface charging process of two 
parallel planar silica surfaces upon decreasing their sep- 
aration, that is the charge regulation. The dependence 
of a on separation, h, in the presence of 0.2 mM Ca 2+ 
at various pH is shown in Fig. 2] At low pH, a mono- 
tonic decrease of a with decreasing h, in conjunction with 
an increased osmotic repulsion (not shown), is found in 
qualitative agreement with previous mean-field studies 
[32J, [33J. This result can be rationalized as a balance 
between entropy and energy. That is, when pH is well 
below p-Koj ol is small and the free energy is dominated 
by the entropic term causing the particles to repel. As a 
result, the system minimizes its free energy by a reduc- 
tion of the surface ionization. Conversely, atypical a vs 
h curves are predicted for surfaces immersed in divalent 
electrolytes at large pH values (pH > ~pKo), where a is 
found to increase when shortening h. The degree of dis- 
sociation reaches a maximum at h max w 6 A, whereafter 
it decreases and annihilates at contact. The maximum in 
a corresponds exactly to the minimum of the potential 
of mean force (not shown). In this case, the energy (in- 
duced by the ion-ion correlations) dominates, which leads 
to an attraction between the particles and, subsequently, 
to an increase of the surface ionization for short h, which 
in turn strengthen the attraction. Close to contact the 
entropy regains in importance and eventually causes a 
to collapse. However, when the Ca 2+ concentration is 
increased to more than 10 mM (not shown), the rise in 
a is dramatically strengthen while the drop at close to 
contact is hardly measurable. In other words, the sur- 
faces retain and even increase their charges all the way 
to contact. 
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This charge regulation behaviour may partly explain 
why Ca 2+ promotes the experimentally observed disso- 
lution of amorphous and crystalline silica surface under 
stress ■ Indeed, as recently rationalized by Dove et al. 
Q , the dissolution proceeds through a nucleation process 
following the same mechanistic theory as developed for 
growth and it is shown to be promoted by surface ioniza- 
tion. 

Ion-ion correlation also seems to be relevant for mem- 
brane selectivity. As an example, we have calculated the 
L-type Ca-channel selectivity at various pH in a reservoir 
containing always 100 mM Na + to which is added either 
1 mM or 100 fiM Ca 2+ . The Ca-channel was modeled as 
an infinite slit pore with h — 9 A, with walls decorated 
with titratable carboxylic groups, p-Ko = 4.8 and p s = 
3/nm 2 . The selectivity, £, was calculated as the ratio 
between the average calcium and sodium concentrations. 
At pH 7 and 100 pM of Ca 2 + we found £=0.96 and a = 
0.70, while at 1 mM Ca 2+ the values increased to £=3.7 
and a — 0.84. Boda et al. [l2| tried to mimick the selec- 
tivity with a constant ionization of the -COOH groups, 
but were then forced to decrease the relative dielectric 
permittivity of the protein to 10 in order to fit the ex- 
perimental data. When pH is reduced to 5, £ drops to 
0.10 and 0.15 at low and high calcium content, respec- 
tively in agreement with the observed loss in membrane 
selectivity at low pH [3s| . 

To summarize, we have shown by confronting Monte 
Carlo simulations with independent experiments that the 
dominating interaction which controls the charging pro- 
cess and charge reversion of silica in multivalent elec- 
trolyte is of purely electrostatic origin and strongly de- 
pendent on ion-ion correlations. In particular, we have 
demonstrated that the charging process and charge re- 
versal are intimately related. We have also demonstrated 
that the accepted view on charge regulation as a mono- 
tonic drop in \a\ with decreasing separation is qualita- 
tively wrong in highly coupled systems. Instead, |er| is 
found to increase with decreasing h. This has profound 
influence for the stability of colloidal particles and for the 
early setting of normal Portland cement. 

We thank R. Kjellander and S. G. Lemay for helpful 
discussions. M. Dove and F.H.J, van der Heyden are 
greatly acknowledge for the kind permission to use their 
experimental data. 
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